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35392 GieBen, Germany 
Using the cluster expansions for n-point Green functions we derive a closed 
set of dynamical equations of motion for connected equal-time Green functions 
by neglecting all connected functions higher than 4*'' order for the A"l> ''-theory in 
■ 1 + 1 dimensions. We apply the equations to the investigation of spontaneous 

ground state symmetry breaking, i.e. to the evaluation of the effective potential 

0\ ■ 

at temperature T — 0. Within our momentum space discretization we obtain a 
CL(' second order phase transition (in agreement with the Simon-Griffith theorem) and 

^' ,2 

X^/y ' critical coupling of Acrit/4m — 2.446 as compared to a first order phase transition 

' and Xcrit/irn^ ~ 2.568 from the Gaussian effective potential approach. 



I. INTRODUCTION 



> 
in 

^T) • In the last years there have been many attempts to develop non-perturbative methods for quan- 



tum field theoretical (QFT) problems such as the infrared behaviour of QGD. Many of the non- 



perturbative methods used so far, such as the coupled cluster expansion originate from 



standard many-body theory and use states in the Fock representation, thus they have the general 



Q^, problem of Fock representations being disjunct to the Hilbert space of QFT |^,^. 

In this work we also adopt a successful method from standard many-body theory, which how- 
• ever is formulated in terms of Green functions and therefore of a genuine field-theoretical nature. 
^ ' We propose an approach along the line of n-body correlation dynamics which describes the 
' propagation of the system in terms of equal-time Green functions (i.e. we use the density matrix 
formalism). The basic strategy of correlation dynamics is the truncation of the infinite hierarchy 
of equations of motion for n-point Green functions via the use of cluster expansions, i.e. the expan- 
sion of Green functions in terms of connected Green functions (correlation functions). In this work 
we go up to the connected 4-point level, i.e. we include the 2-, 3- and 4-field correlation functions. 
For relativistic systems with a nonlocal interaction our method does not guarantee a covariant de- 
scription, since retardation effects can not be included due to the equal-time formalism. There is, 
however, no objection against the use of an equal-time formalism in the case of a local relativistic 
field theory as long as one consistently considers Green functions containing canonically conjugate 
field momenta as well as the fields themselves. 
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This paper is organized as follows: In section || we describe the derivation of the correlation 
dynamical equations of motion for $*-theory in 1 + 1 dimensions; the final set of equations itself is 
shifted to appendix ^ in view of its length. Section III is devoted to the application of the method 
to the evaluation of the effective potential at zero temperature within various limits and thus to 
the investigation of ground state symmetry breaking. 



II. CORRELATION DYNAMICS FOR $l+i-THEORY 



We consider the Lagrangian density 
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in one space and one time dimension, which corresponds to the Hamiltonian 

"2 



i7 = i \ dx 
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(11.1) 



(11.2) 



where t: — dt(j) and mo is the bare mass of the scalar field 0. In order to evaluate the effective 
potential of the theory at zero temperature, i.e. the minimum of the energy density for a given 
magnetization (</>), we decompose into a classical and a quantum part according to 



where $o is a real constant and assume 
($)=0, i.e. (0)=$o- 

Since $o is a constant, we have 
t: = dt(t) = dt<i> = n . 



(11.3) 



(II.4) 



(11.5) 



The Hamiltonian, expressed in terms of the classical part $o and the quantum part <& of (p, then 
reads 



H^^ I dx 



+ ( II ) + (2$o"i§ + 2X^1) $ + (m^ + 3A$2) $2 ^ 2A$o*^ 



-iA$4 + ^2^2 + iA$^' 



(11.6) 



For the time evolution of <i> and 11 we obtain with (11.6) and the canonical equal-time commutation 
relations 



t), $(y, t)] [Il{x, i), n(y, t)] = , Mx, t), n(y, t)] = i6{x - y) 
by means of the Heisenberg equation: 



(II.7) 



dtH = - {^oml + A$g) + {d^ -ml- 3A$g) $ - 3A$o*^ - 



(II., 
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In analogy to ([1.8) we get the equations of motion for equal-time operator products of $ and H, 
e.g. 

dt , i)$(a;2 , i)) = n(xi , t)$(x2 , i) + , t)U{x2 , t) , 

dt {n{xi,tMx2,t)) = [- ($0^2 + A$ij) + {dl ~ml^ 3A$g) t) 

-3A$o$^(a:i,t) - A$3(a;i,t)] $(x2,t) + n(a;i,t)n(x2,0 , (11.9) 

After taking their expectation values, the equations of motion for all possible n-point equal-time 
products of $ and 11 comprise an infinite hierarchy of equations of motion for equal-time Green 
functions. The analogon of this in nonrelativistic many-body theory is the BBGKY|^ density matrix 
hierarchy . Since we are considering a local field theory, the equations of motion for the Green 
functions contain no retardation integrals, which implies that working in an equal-time limit is 
sufhcient for describing the propagation of the system as long as one considers all Green functions 
containing the fields as well as their conjugate momenta. 

For practical purposes, the infinite hierarchy of coupled differential equations of first order in 
time has to be truncated. This is done using the cluster expansions for n-point Green functions [Q, 
i.e. their decomposition into sums of products of connected Green functions, which works as long 
as the system is in a pure phase [^,^. The explicit form of the cluster expansions can be derived 
from the generating functionals of full and connected Green functions, Z[J^a] and W^[J, cr], given 
by 

Z[J,a] = TrjpT [e*/'''*(^(*)*("^+"(*)n("»] } and Z[J,a] = e^^''^"^ , (11.10) 

respectively, where T is the time ordering operator, p is the statistical density operator describing 
the pure or mixed state of the system (Trp = 1) and x = [x, t). We start with the cluster expansions 
for the time-ordered Green functions with different time arguments: 

{Hx)) = ($(£)), , (n(x)) = (n(x)), , (11.11) 

(T$(xi)$(x2)) = lim -J—-J—e'^y-J^-\ 

Jm^O iSJ(Xi) lSJ[X2) 

= lim -^^^\(-^W[J,a])e'^i''^ 

J.a^O l5J{xi) [V*'^'^(^2) V 



6 S 
j,a^a 1^ yidJ{xi) i5J{x2) 



= lim 



= (T$(xi)<i>(i2)), + (<i>(ii))(<i>(£2)) , (11.12) 

where {■)c denotes the connected part of the expectation value. Analogously we obtain 

(rn(a;i)$(i2)) = (rn(£i)$(£2))c + (n(.Ti))($(x2)) , 



^Born Bogoliubov Green Kirkwood Yvon 
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(Tn(fi)n(i2)) = (Tn(fi)n(i2))c + (n(fi))(n(i2)) , 



(r$(£i)$(£2)$(x3)) = (T$(X1)$(X2)$(£3)), + (T$(X1)$(X2)),($(£3)) 

+(r$(£i)$(£3)),($(£2)) + (r$(£2)$(£3)),(*(£i)) + ($(^i))($(:£2))($(i3)) , ... . 

(11.13) 

The expressions for equal-time Green functions are obtained by taking the weU-defined equal-time 
limit which yields the appropriate operator ordering in the cluster expansions. We arrive at 

{<f{x)) ^ mx))^ , {u{x)) = (n(x)), , 



($(X1)$(X2)) = ($(X1)$(X2)), + ($(X1))($(X2)) , 

(n(xi)$(a;2)) = (n(xi)$(a;2))c + (n(xi))($(x2)) , 
(n(xi)n(a;2)) - (n{xi)U{x2))c + (n(a;i)}(n(x2)) , 



($(xi)$(a;2)$(x3)) - (*(a;i)$(x2)$(a;3)),-F($(xi)$(a;2)),($(a;3)) 

+ {Mxi)^Xs)),mX2)) + mX2Mxs)),mX,)) + ($(X1))($(X2))($(X3)) 



(11.14) 



where all (equal) time arguments have been suppressed. In view of their length the cluster expan- 
sions for the other Green functions required for our calculations are not explicitly given here, but 
e.g. can be found in |l0| . 



Due to equations ( II. 4 ) and ( 11.5 ) the 1-point functions in all cluster expansions are now assumed 
to vanish, which in the end leads to the same result as considering the cluster expansions for Green 
fmictions containing the original fields and tt and setting (</)) = $o £ind (vr) = 0. 

The cluster expansions then are truncated by neglecting all connected n-point Green functions 
with n > N; in our case with n > 4. Inserting the truncated cluster expansions into the equations 



of motion of type (II.9) up to the equations for the 4-point functions leads to a closed system 
of coupled nonlinear equations for the connected Green functions (i.e. a system of correlation 
dynamical equations) up to the 4-point level in analogy to the n-body correlation dynamics in 
nonrelativistic many-body theory For this straightforward but somewhat tedious derivation, 
which in view of its length is not explicitly given here, all truncated cluster expansions up to 
the 6-point level are required, since the highest order Green functions appearing in the hierarchy 



equations (II. 9) up to the 4-point level are the 6-point functions. 

The resulting equations of motion for the connected Green functions still have to be renormalized. 
<l?^-theory in 1 -|- 1 dimensions is superrenormalizable and only requires a mass renormalization. 
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There is only one mass counterterm due to the divergent tadpole diagram comprising the contri- 
bution to the selfenergy of lowest order in the coupling constant. This mass counterterm can be 
evaluated analytically either in the framework of perturbation theory or, equivalently, by normal 
ordering the Hamiltonian with respect to the perturbative vacuum^. We obtain pi] ]: 

ml = ir? + 8jr? , S-m? ^ -3AA(0) (11.15) 

with 



dp e 



The logarithmically divergent mass counterterm ( 11.15 ) can be analytically removed from the equa- 



tions of motion by normal ordering all operator products with respect to the perturbative vacuum 
and thereby splitting off the short-distance singularities of the free equal-time Green functions. 
The normal ordering within the connected Green functions only takes place within the 2-point 
functions, since in the other connected n-point functions all field operators commute with each 
other. For the 2-point functions we have: 

($(xi)$(a;2))c = (: «'(xi)$(x2) :)c + A(a;i - X2) , 
(n(xi)n(a;2)), = (: n(xi)n(a;2) :)c - {d^^ - m^) A(xi - x^) , 
(n(xi)$(a;2))c = (: n(xi)$(a;2) :)c - ^<5(xi - X2) , 

mxi)U{x2))c = {■■ *(xi)n(a:2) :)c + ^<5(xi - x^) . (11.17) 



After inserting ( [1.15 ) and ( II.17| ) into the equations of motion for the connected Green functions. 



all terms containing a factor A(0) mutually cancel out, leading to a closed set of renormalized 
equations for the normal ordered connected Green functions. 

The final step is to transform the renormalized correlation dynamical equations of motion from 
coordinate space to an arbitrary single particle basis in order to simplify their numerical integration. 
In this respect we expand the field $ and its conjugate momentum 11 according to 

$(x) = ^*aV'a(a;) , n(a:) = J^'^c.tpaix) with / dx C(a;)V'/3(a;) = S^p . (11.18) 

a a 

For the corresponding equal-time Green functions wc then have 

mxi)^{x2))c = '^{^a^l3)c'4'a{xi)ij0{x2) = (: ^{xi)^X2) :) c + ^{xi - X2) 
a(3 

= $Q*/3 ■)c1pa{xi)iJi3(x2) + '^^at3lpa{xi)tpi3{x2) , 

a/3 a/3 



^In addition to introducing the mass counterterm, one of course has to subtract the constant infinite 
zero-point contribution from the Hamiltonian, which however does not enter the Heisenberg equations. 



5 



a/37 

($(xi)$(a;2)«>(x3))c = ^($a$/3$7)cV'a(a;i)V'/3(a;2)V'7(a;3) , - ■ (11.19) 

a/37 

The corresponding expressions for the other Green functions are obtained analogously. By inserting 
the expansions ( II.19D into the renormalized equations of motion and projecting out the matrix 



elements with respect to the single particle basis we obtain the final result for the correlation 
dynamical equations of motion for the $1^]^ -theory, which are denoted by $*CI? Correlation 
Dynamics) furtheron, and can directly be used for a numerical integration. In view of their length 
the $*CZ3 equations are shifted to appendix 0. 

III. APPLICATION TO GROUND STATE SYMMETRY BREAKING 

In this section we apply the <I>^CZ? equations to the determination of the effective potential of 
the $^^j^-theory at zero temperature^ and thereby investigate the spontaneous breakdown of the 
symmetry under the discrete transformation —>■ —0 for values of the coupling exceeding a critical 
value, which manifests itself in a nonzero ground state magnetization (</>). 

Since we also want to study the influence of the different connected n-point functions, we in- 
troduce 4 limiting cases of the correlation dynamical equations; these are denoted as $'*CD(2), 
^^CD{2,5), <i>''C£>(2,4) and ^^CD{2,S,A), where the numbers in parantheses are the orders of 
connected n-point fmictions that are taken into account (i.e., for instance in the ^"^00(2, 4) case 
the connected 3-point functions are set equal to zero and $4CD(2, 3, 4) denotes the original <I>^CD 
approximation discussed in section ||) . 



In order to integrate equations ([A.l|) - (A. 12) numerically, we choose plane waves in a one- 



dimensional box with periodic boundary conditions as a single particle basis, i.e. we work in 
discrctized momentum space. In order to select an appropriate box size for a given renormalized 
mass, we compare the numerically obtained GEPQ solution for the ground state configuration 
within the discrctized system to the analytically accessible GEP solution in the continuum limit 
(for the GEP approximation, see appendix |b| or ]ri|-]l3[). Due to the large amount of computer 
time needed for ^"^CD calculations we have to compromise between a good momentum space 
resolution and the convergence with a minimum number of plane waves. It has proven to be most 
effective to choose a box size of 100 fm for a renormalized mass of 10 MeV. In fig. |l] we show the 
ground state magnetization (0) = $0 a-s a function of the dimensionless coupling X/Am^ for these 
parameters, where different numbers of single particle states have been taken into account. The 
vertical line shows the position of the critical coupling for symmetry breaking in the continuum 



^i.e. the ground state energy density in a subspace with a given magnetization 
^Gaussian Effective Potential 
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limit. For practical purposes we will always use the 15 lowest lying plane waves in the following 
calculations, unless explicitly stated otherwise. Within the GEP approximation we then obtain a 
critical coupling of Acrit/4m^ = 2.568 for the discretized system as compared to a critical coupling 
of Xcrit/^'m? — 2.5527045 in the continumn limit. 

Since the ^'^CD equations only describe the propagation of equal-time Green functions for a 
given initial configuration, they cannot directly be applied to the evaluation of static equilibrium 
properties of the system. As in the case of correlation dynamics for nonrelativistic many-body 
theory there is no easy access to the stationary solutions of the equations, and moreover it is 
not clear, in how far additional constraints have to be imposed on the subspace of stationary 
configurations in order to select only the physical solutions. 

We therefore use a different approach for the evaluation of equal-time ground state Green func- 
tions within the $^CZ3 approximation. Starting with the trivial exact ground state configuration 
for A = and a given fixed value of $0 as an initial condition, we continuously switch on the 
coupling while propagating the system in time. The time-dependence of the coupling is chosen to 
be linear with \ ~ [3t and (3 = const. 

In fig. ^ the energy density obtained within this time-dependent method is shown as a function 
of the coupling for different values of (3 for $0 — 0.4 and all 4 limiting cases of the correlation 
dynamical equations. In all cases an asymptotic curve is approached with decreasing (3, i.e. when 
the coupling is switched on more slowly. This indicates that in the limit /3 ^ the whole process 
becomes fully adiabatic, i.e. the system will time-dependently follow the trajectory of the ground 
state as a function of the coupling. For the <I>^C£'(2) approximation, taking into account only the 
connected 2-point functions, which is the field theoretical analogon to time-dependent Hartree-Fock 
theory, we have direct access to the static ground state solutions since its stationary limit is simply 
given by the GEP approximation. The corresponding curves are shown in the upper left part of fig. 

the GEP is displayed as the lower one of the two solid lines. Indeed, the asymptotic curve of the 
time-dependent method is identical to the GEP solution. In addition, we have checked numerically 
that for all 4 limiting cases of correlation dynamics the system increasingly equilibrates along its 
trajectory with decreasing (3, i.e. if we stop increasing A at some point in time, the system will 
remain in its present state when propagated further. 

In choosing a finite value for (3 we have to compromise between a good convergence of our time- 
dependent method and a minimum of computer time we want to invest; the curves in the two 
following pictures have therefore been evaluated with /3 — 5MeV^c/fm. 

The ground state energy density for various values of $0 is plotted versus the coupling A/4m^ 
in fig. H for all 4 limiting cases of correlation dynamics, where for the reasons mentioned above 
the <I>^CZ?(2) approximation has been replaced by the GEP approximation. The GEP and the 
<I>^C£'(2,4) approximation each predict a first order phase transition, since the first curve to 
intersect the $0 = (i.e. the symmetric phase) energy density in the GEP case has $0 = 0.6 and 
in the $''C£'(2,4) case has $0 = 0.8; i.e. they both lead to a finite value of $0i implying that 
the vacuum magnetization has to jump to that value discontinuously at the critical value of the 
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coupling constant (compare fig. |i|for the GEP). In contrast to that, the ^^CD{2, 3) approximation 
and the <i>'*C£'(2, 3, 4) approximation each predict a second order phase transition, since in both 
cases all curves with $o 7^ intersect the curve with $o = in the correct order, thus enabling 
the vacuum magnetization to increase continuously once the critical coupling is exceeded. 

At this point it is useful to recall, that there is a rigorous mathematical proof of the statement 
that there can be no first order phase transition in the $|^j^-theory (in the continuum limit) Jlif . 
This proof is based on the Simon-Griffith theorem |]l5| , which in turn is obtained by considering 
the field theory as a proper limit of a generalized Ising model. 

Thus we conclude, that in order to describe the order of the phase transition correctly, the 
inclusion of the connected 3-point function is required. A simple "geometrical" explanation of this 
fact will be given in the discussion below. 

In fig. ^ the effective potential is plotted as a function of the magnetization $o for various values 
of the coupling. The same data have been used as for the previous figure. In this representation, 
the first order nature of the phase transition in the GEP and the $'*CZ3(2, 4) case can be seen from 
the fact, that for these approximations there is a maximum in the effective potential between the 
minimum at $o = and the second minimum on the right hand side (for the GEP this minimum 
is not very pronounced); thus at the critical coupling the second minimum has to be located at 
a finite value of $o- In the case of the $^C£'(2,3) and the $'*CZ3(2, 3, 4) approximation there is 
no such intermediate maximum, and as the critical coupling is approached from above, the second 
minimum is shifted to the left towards $0 = 0. 

For $0 ^ c» the system will approach the classical limit, i.e. with decreasing $o higher order 
correlations increasingly become important; the inclusion of correlations always lowers the energy 
density, since a larger configuration space is opened up for the system. However, for obvious 
mathematical reasons the inclusion of the 2-point function alone cannot lower the energy density 
at $0 = 0; unlike the other n-point functions of even order the 2-point function therefore is most 
important in the region on the right hand side in the plots of fig. ^. The 4-point function, as 
the highest order correlation included in this work, is most important around $o = 0, i.e. for 
small classical field strengths. The intermediate maximum now appears in between the domains 
governed by the 2-point function on the right and the 4-point function on the left, an inclusion of 
higher correlations of even order would most probably lower the energy density in a region even 
more concentrated around $o = and therefore not cure this problem. The inclusion of connected 
n-point functions of uneven order cures the problem, since for symmetry reasons these have to 
vanish as $o ^ 0. The 3-point function therefore obviously assumes its main importance exactly 
in the region in between the domains of the 2- and the 4-point function, it "bends down" the 
intermediate maximum. 

Finally, in table | we give the values for the critical coupling extracted from our calculations for 
different numerical parameters and the 4 different limiting cases of $*C-D. The critical couplings 
Acrjt/4?n^ obtained by other authors, that are also using non-perturbative techniques adopted from 
standard many-body theory, are as follows: Jl6| , p7t obtain 1.829 or 1.375, respectively, using the 
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discretized light-front quantization method, [^8[ obtain 1.72 using a varitional approximation with 
trial states which are quartic rather than gaussians (quadratic exponentials) causing an inclusion 
of the 3-field correlation amplitude in the language of the coupled cluster method 0] , ||] obtain the 
estimate 0.95 < A/4m^ < 2.15 using an improved version of the coupled cluster expansion method, 



and using second order perturbation theory in the residual interaction leads to 1.14 |19|. All of 
these authors obtain a second order phase transition in agreement with the Simon-Griflith-theorem. 



approximation 


No. of states 


/3 [MeV^c/fm] 


Acrit/4m^ 


order of phasetransition 


GEP 


15 




2.568 


first 


$''CLi(2,3) 


11 


5.0 


1.679 


second 


13 


5.0 


1.649 


15 


5.0 


1.629 


17 


5.0 


1.613 


19 


5.0 


1.601 


$''CD(2,4) 


15 


5.0 


3.81 


first 




15 


5.0 


2.446 


second 



TABLE I. Values for the critical coupling for different approximations and numerical parameters. 
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IV. SUMMARY AND OUTLOOK 



This work comprises the first apphcation of correlation dynamics for equal-time Green functions 
to a field theoretical problem, i.e. to the determination of the effective potential in $f^]^-theory. 
After giving a derivation of the corresponding equations of motion, we showed that we are able to 
evaluate equal-time quantities in the interacting ground state for a given vacuum magnetization 
by time-dependently increasing the coupling in an adiabatic process. 

Our numerical results predict a second order phase transition in agreement with the Simon- 
Griffith theorem, as soon as the connected 3-point function is included. However, since the con- 
nected 2-point function and the connected 3-point function alone are not able to lower the energy 
of the symmetric phase, the critical coupling obtained within the ^^CD{2, 3) approximation is too 
low. Going one step further in our expansion, we find that the connected 4-point function in the 
^^CD{2, 3, 4) method increases the coupling to a value, which is only slightly lower than the value 
obtained by the GEP approximation, where however the shape of the effective potential changes 
completely as compared to the GEP. 

In general, the results of this work demonstrate the applicability of correlation dynamics to the 
description of low-energy (ground state) phenomena in local field theories. Since in principle the 
equations are designed to describe the propagation of the system in time with arbitrary initial 
conditions, our method is also a potentially powerful tool for the investigation of non-equilibrium 
properties of relativistic field theories, e.g. the response to external perturbations. 

The present study can be viewed as a first step towards the application of correlation dynamics 
to SU(N) gauge theories, aiming at a non-perturbative description of the infrared behaviour of 
QCD; we already presented the corresponding equations of motion in |ic| ]. 

APPENDIX A: EQUATIONS OF MOTION OF CORRELATION DYNAMICS 

In this appendix we present the renormalized correlation dynamical equations for the normal 
ordered connected Green functions up to the 4-point level, formulated with respect to an arbitrary 
single particle basis set (see section ||). In order to compactify the equations we introduce the 
following abbreviations: 
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t^p = Jdx raix) [dl - - 3A$^) i)0{x) , 
Ua0 = -3A 5^(a|AiA2/3)(: $Ai$a. :)c • 

Ai A2 

The permutation operator interchanging the indices a and (3 is denoted by Pap- Although normal 
ordering only afFccts the connected 2-point functions, wc write out the normal ordering operation 
: • : in all connected Green functions in order to have a uniform notation. The equations then read: 

^^^0 :)e = (: U^^0 + (: $„n^ , (A.l) 



nc^0 :)c = (: H^H^ :)c 

+ 5I(*«A + C^aA)(: ^A$/3 :)c + $^(f/aA - 3A$^5„a)Aa/3 

A A 

-A ^ (a|AiA2A3)(: *Ai*A2*A3*/3 :)c 

Ai A2 A3 

-3A$o Yl ("l^i^2)(: $Ai$A2*/3 :)c , (A.2) 

Ai A2 



n„n^ :)e = (1 + Pa^) ^(iaA + Uc,x){: ^xll0 :)c 
-A(l + P«/3) ^ (a|AiA2A3)(: $Ai$A2$A3n/3 :)c 

Ai A2 A3 

-3A$o(l + Va0) J2 ("l^i^2)(: $Ai$A.n;3 , (A.3) 

Ai A2 



$a$/3*7 :>c = (: n«$/3$^ :)c + (: ^alif3% ■)c + {■■ '^a^ff^j :)c , (A.4) 



nc,$;3$7 :)c = (: Tla^Up^^ :)e + (: Happily :)c 
A 

-3A J2 (a|AiA2A3>(l+P/37)(:$Ai$A2^/3:)c((:$A3^7:)c + AA37) 

Ai A2 A3 

-3A$o 5I("|AlA2){(: $Ai$A2^/3$7 Oc 
Ai A2 

+2((: $A,$/3 :)c + Aa,/5)((: $a.$7 Oc + Aa,7)} , (A.5) 
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+ (1 + ^(iaA + Uax){-- ^aH^^t 

A 

-3A(1 + Vc,0) (a|AiA2A3){(: ^Ai^A^H/s :)c((: ^x,% :)c + AA37) 

Ai A2 A3 

+ (: $Ai$A.$7 :)e(: '^'x^^fi Oc} 

-3A$o(l+7'a/3) ^(a|AiA2){(: ^Ai^a.H^s^^ :)e 

Ai A2 

+2(: $;,,n^ :),((: ^x.^j :)c + Aa,^)} , (A.6) 



-(: n^n^n^ ■.)^ = (l + Vaf3 + Va^) Y,{taX + f/aA)(: ^X^fiHj ■)c 

-3A(1 + P„;3+P«7) I] (a|AiA2A3>(l+P/3^)(:$Ai$A.n^:>c(:$A3n^:)c 

Ai A2A3 

-3A$o(l + Pa/3 + T'aT) (a|AiA2){(: $Ai*A.n^n^ :)e + 2(: $Ain;3 $A.n^ :>c} 

Ai A2 

+ ^A$o(a/37|l) , (A.7) 



dt 



+ (: <i>aUi3%^s ■)c + {■■ 4'a$/3n^*5 :)c + (: $a*/3*7n5 ■)c , 



(A.8) 



A 

-3A ^ (a|AiA2A3){2((: $Ai*/3 :>c + ^X,0){{- ^X^^j ■)c + ^X^jXi- ^X^^S ■)c + ^X^s) 

Ai A2 A3 

+ + $Ai$;3 :)c + Aai/3)(: $a.$A3*7*5 :>c 

+ (: $Ai$A2*/3 :)c(: $A3^7^5 Oc]} 

-6Ai>o ^ (a|AiA2)(l + ^^37 + P^sM *Ai*/3 Oc + $a2^7^^ •)c > (A.9) 

Ai A2 
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+ (1 + Va0) Y^ito^X + Uax){-- ^X^0^j^S ■)c 

X 

-3A(1 + Pa0) (a|AiA2A3){2(: ^^.U^ :),{{: $^.$7 ■}c + ^x.-yM <^Xs<^S :)c + ^Xss) 

Ai A2 A3 

+(: ^x.'ilff :)c(: $a,$A3^7^5 Oc + (1 + V^sXi'- $Ai$7 :)c + Aa,^)(: ^x^^Xs^/S^S :)c 
+(1 + Vp^ + Vps){: ^Ai^A^n/j :)c(: ^X,^^^5 Oc} 

-6A$o(l+7'„^) X^(a|AiA2){(: ^a^H^ :)c(: *a.*7*5 

Ai A2 

+ (1 + V^5){{: $Ai$7 Oc + Aai7)(: $A2n/3$5 :)c} , (A.IO) 



-(: n^n^n^*^ = (: n^n^n^n^ 

+ (l + ^a/3+^a7)13(iaA + C/aA)(: '^x'^p^^'^S :)c 
A 

-3A(1 + + Va^) J2 (a|AiA2A3){2(: ^a^H^ :)c(: ^X.^^ :>c((: ^X,^5 :)c + AA35) 

Ai A2 A3 

+(1 + Vi3j){: ^Ai^A^n/s^-a :)c(: $A3n7 :)c + (: ^Ai^A^H^H^ :)c((: $A3^5 :)c + Aa35) 
+{1 + V0^ + V05){- ^Ai^A^n/j :)c(: ^x.U^'^S :)c} 

-6A$o(l + ?'a/3 + 7'a7) I^(a|Al^2){(l+P/37)(: ^Xi'^fi OcO $A2n7$5 :)c 
Ai A2 

+(: ^AiH^n^ :)c((: ^A^^i Oc + Aa,^)} 

+ ^A^(a/37|A)((: ^A'i?^ :)c + Aa^) , (A.ll) 

A 

n^n^n^n^ = (i + + + p^s) X^(iaA + Uax){- ^An^n^n^ 

-3Xil + V„0 + Vaj+Vas) J2 {a\XM3){2{:^x,Tli3:)c{:^x,Tl^:)c{:^Xs^s:)c 

Ai A2 A3 

+ + $Ain/3 *A.*A3n7n5 :)c + (: $Ai$A.n^ :)c(: $A3n7n5 :)c]} 

-6A$o(l + ^'a/J + Vaj + Vas) ^ (a|AiA2>(l + P/37 + 7'/35)(: ^A^n^ :)c(: ^X.^^^S :)c 

Ai A2 

+ ^A(l + Va0 + Va-y + Vas) ^(/37'5|A)(: ^X^a :)c • (A.12) 

A 
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APPENDIX B: THE GAUSSIAN EFFECTIVE POTENTIAL FOR $^+i-THEORY 



The Hamiltonian for the <i>f^]^-system is given by 



H= - I dx 



with rriQ = m? + 6m?, 5m? = — 3AA(0) (see section Let furthermore 

= (27r)"5 j dp {2ujp)~'^ + aLp 

ap\0) = 0, 



$0 + ^', 



where |0) is the perturbative vacuum. As in section |l|, $0 denotes the constant vacuum magneti- 
zation, i.e. the ground state expectation value of (j). 
The GEP approximation consists in the ansatz 



Hf) = - I dx 



(B.2) 



for the variational Hamiltonian, i.e. the interacting system is approximated by a free system with 
an effective mass M, which serves as a variational parameter. One then has to minimize the 
expectation value of H with respect to the ground state of Hq. 

This method is equivalent to a BCS calculation, where the variational wavefunction is given by 
a boson pair condensate 

|*> =Mexp{-S)\0) , S^^Jdp fi{p)alalp ; (B.3) 

here Af denotes a normalization constant and fi{p) is the variational parameter. |^) is a quasipar- 
ticle vacuum, i.e. fop|^) — 0, where the bp can be obtained from the Op by means of a Bogoliubov- 
transformation. We have 



= (27r)"2 J dp {2np)-^ 

with np = ujp^±^. 

After minimizing the energy functional, the GEP ansatz yields the equation 

dp 1 1 



(B.4) 



3A $^ 



27r V 2Q,p 2ujp 



(B.5) 



with Vlp ~ ^/p^ + for the effective mass M. Equation ( |B.5| ) is identical to the HFDS0 equation 
and the Gap equation following from the BCS ansatz [0. 



For known effective mass M, the equal-time Green functions can be evaluated via (B.4) as 
expectation values with respect to the quasiparticle vacuum The value for the energy den- 
sity functional as a function of $0 (obtained with the method described above) is the Gaussian 
approximation for the effective potential of the theory. 



^Hartree Fock Dyson Schwinger 
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Figure Captions 

fig.l: Vacuum magnetization as a function of the coupling in the GEP approximation for 
different numbers of basis states. 

fig. 2: Ground state energy density during the time-dependent process as a function of the 
coupling for a given vacuum magnetization of $o = 0.4 and different values of /3; 
upper left: <i>'^CD{2), upper right: ^'^CD{2,3), lower left: <i>^CD{2,4:), lower right: 
$'*C£'(2, 3, 4); in addition the static GEP curve is displayed in the upper left plot. 

fig. 3: Ground state energy density as a function of the coupling for diflFerent vacuum magne- 
tizations; all curves except the one in GEP approximation have been evaluated time- 
dependently with f3 = 5 MeV^c/fm; upper left: GEP, upper right: $^CD(2,3), lower 
left: $4C£)(2, 4), lower right: $'^CL>(2, 3, 4). 

fig. 4: Effective potential as a function of the vacuum magnetization for different values of 
the coupling; all curves except the one in GEP approximation have been evaluated 
time-dependently with /? = 5 MeV^c/fm; upper left: GEP, upper right: $*C£)(2,3), 
lower left: $^C£>(2,4), lower right: $^C£)(2, 3, 4). 
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